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Abstract 

Solar radiation models for predicting the average daily and hourly global radiation, beam 
radiation and diffuse radiation are reviewed in this paper. Seven models using the Angstrom- 
Prescott equation to predict the average daily global radiation with hours of sunshine are 
considered. The average daily global radiation for Hong Kong (22.3°N latitude, 114.3°E 
longitude) is predicted. Estimations of monthly average hourly global radiation are discussed. 
Two parametric models are reviewed and used to predict the hourly irradiance of Hong Kong. 
Comparisons among model predictions with measured data are made. 0 2001 Elsevier 
Science Ltd. All rights reserved. 


1. Introduction 

A knowledge of the local solar-radiation is essential for the proper design of 
building energy systems, solar energy systems and a good evaluation of thermal 
environment within buildings [1-6]. The best database would be the long-term 
measured data at the site of the proposed solar system. However, the limited cover¬ 
age of radiation measuring networks dictates the need for developing solar radiation 
models. 

Since the beam component (direct irradiance) is important in designing systems 
employing solar energy, such as high-temperature heat engines and high-intensity 
solar cells, emphasis is often put on modelling the beam component. There are two 
categories of solar radiation models, available in the literature, that predict the beam 
component or sky component based on other more readily measured quantities: 
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Nomenclature 

A 

Apparent extraterrestrial irradiance (W m -2 ) 

a\, 02, a-i 

Coefficients defined in Eqs. (33), (50) and (68), respectively 

B 

Overall broadband value of the atmospheric attenuation- 

b u b 2 , b 3 

coefficient for the basic atmosphere (dimensionless) 
Coefficients defined in Eqs. (33), (50) and (68), respectively 

C 

Ratio of the diffuse radiation falling on a horizontal sur- 

c„ 

face under a cloudless sky 4 to the direct normal irradia¬ 
tion at the Earth’s surface on a clear day 4 (dimensionless) 
Clearness number, the ratio of the direct normal irradiance 

C 1 ,C2,C 3 ,C' 3 ,C 4 ,C' 4 

calculated with local mean clear-day water vapour divided 
by the direct normal irradiance calculated with water 
vapour according to the basic atmosphere (dimensionless) 
Coefficients defined in Eqs. (81)—(92) 

A 

Aerosols scattered diffuse irradiance after the first pass 

D m 

through the atmosphere (W m -2 ) 

Broadband diffuse irradiance on the ground due to the 

A 

multiple reflections between the Earth’s surface and its 
atmosphere (W m -2 ) 

Rayleigh scattered diffuse irradiance after the first pass 

d\, d 2 , d 3 , dn, d$, df, 

through the atmosphere (W m -2 ) 

Coefficients defined in Eqs. (81)—(92) 

E 0 

Eccentricity correction factor of the Earth’s orbit (dimen- 

F c 

sionless) 

Fraction of forward scattering to total scattering (dimen- 


sionless) 

Angle factor between the surface and the Earth (dimen- 

Fss 

sionless) 

Angle factor between the tilt surface and the sky (dimen- 

G b 

sionless) 

Monthly average daily beam irradiance on a horizontal 

G c 

surface (MJ m -2 ) 

Monthly average daily clear-sky global radiation incident 

G 0 

on a horizontal surface (MJ m -2 ) 

Monthly average daily extraterrestrial radiation on a hor- 

G t 

izontal surface (MJ m -2 ) 

Monthly average daily global radiation on a horizontal 

h 

surface (MJ m -2 ) 

Elevation above sea level (km) 

4,4 

Total beam irradiance on a horizontal surface (W m -2 ); 

4,4 

hourly beam irradiance on a horizontal surface (MJ m -2 ) 
Total diffuse irradiance on a horizontal surface (W m -2 ); 
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hourly diffuse irradiance on a horizontal surface (MJ 

In, In 

m -2 ) 

Direct normal irradiance (W m -2 ); hourly normal irra¬ 
diance (MJ m -2 ) 

ioJo 

Extraterrestrial radiation (W m -2 ); hourly extraterrestrial 
radiation (MJ m -2 ) 

Ir 

Reflected short-wave radiation (W m -2 ) 

j sc , /sc 

Solar constant, 1367 W m -2 ; solar constant in energy unit 
for 1 h (1367x3.6 kJm- 2 ) 

h, h 

Global solar irradiance (W m -2 ), hourly global solar irra¬ 
diance (MJ m -2 ) 

k b 

Atmospheric transmittance for beam radiation (dimen¬ 
sionless) 

k D 

Atmospheric diffuse coefficient (dimensionless) 

k d 

Diffuse fraction (dimensionless) 

k t 

Clearness index (dimensionless) 

4o 

Aerosol optical-thickness (dimensionless) 

lo z 

Vertical ozone layer thickness (cm) 

ffio 

Rayleigh optical-thickness (dimensionless) 

m a 

Air mass at actual pressure (dimensionless) 

m r 

Air mass at standard pressure of 1013.25 mbar (dimension¬ 
less) 

N 

Day number in the year (No.) 

P 

Local air pressure (mbar) 

Po 

Standard pressure (1013.25 mbar) 

4rat 

Downward radiation from the atmosphere due to long¬ 
wave radiation (W m -2 ) 

RH 

Monthly average daily relative humidity (%) 

r t 

Ratio defined in Eq. (53) (dimensionless) 

S 

Monthly average daily sunshine-duration (h) 

Sf 

Monthly average daily fraction of possible sunshine, 
known as sunshine fraction (dimensionless) 

So 

Maximum possible monthly average daily sunshine-duration 

(h) 

S 0 m 

Modified average daily sunshine duration for solar zenith 
angle 0 Z greater than 85° (h) 

T 

Monthly average daily-temperature (°C) 

/’at 

Atmosphere temperature at ground level (K) 

/dew 

Dew-point temperature (°C) 

r L 

Linke turbidity-factor defined in Eq. (112) 

/sky 

Apparent sky-temperature (K) 

Ut 

Pressure-corrected relative optical-path length of pre- 
cipitable water (cm) 
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u 3 

Ozone relative optical path length under the normal 
temperature and surface pressure NTP (cm) 

Vis 

Horizontal visibility (km) 

w 

Precipitable water-vapour thickness reduced to the stan¬ 
dard pressure of 1013.25 mbar and at the temperature T of 
273K (cm) 

wf 

Precipitable water-vapour thickness under the actual con¬ 
ditions (cm) 

A 

Change of 

r 

Day angle (in radians) 

a 

Solar altitude (degrees) 

Puh 

Angstrom turbidity-parameters (dimensionless) 

S 

Declination, north positive (degrees) 

Sc 

Declination of characteristics days, north positive (degrees) 

£at 

Atmosphere emittance for long-wave radiation (dimen¬ 
sionless) 

<P 

Latitude, north positive (degrees) 

K Z 

Reflectance of the foreground (dimensionless) 

X 

Wavelength (pm) 

e 

Angle of incidence; the angle between normal to the sur¬ 
face and the Sun-Earth line (degrees) 


Zenith angle (degrees) 

Pa 

Albedo of the cloudless sky (dimensionless) 

Pc 

Albedo of cloud (dimensionless) 

Pg 

Albedo of ground (dimensionless) 

a 

Stefan-Boltzmann constant (5.67xl0 -8 W m -2 K -4 ) 

T r , T 0 , T g , T w , T a 

Scattering transmittances for Rayleigh, ozone, gas, water 
and aerosols (dimensionless) 

T aa 

Transmittance of direct radiation due to aerosol absorp- 
tance (dimensionless) 

T as 

Transmittance of direct radiation due to aerosol scattering 
(dimensionless) 

00 

Hour angle, solar noon being zero and the morning is 
positive (degrees) 

0)1 

Hour angle at the middle of an hour (degrees) 

0) o 

Single-scattering albedo (dimensionless) 

(Os 

Sunset-hour angle for a horizontal surface (degrees) 

f 

Function 

% 

Tilt angle of a surface measured from the horizontal 
(degrees) 

< > 

Yearly average 




L.T. Wong, W.K. Chow / Applied Energy 69 (2001) 191-224 


195 


• Parametric models 

• Decomposition models 

Parametric models require detailed information of atmospheric conditions. 
Meteorological parameters frequently used as predictors include the type, amount, 
and distribution of clouds or other observations, such as the fractional sunshine, 
atmospheric turbidity and precipitable water content [5,7-14], A simpler method 
was adopted by the ASHRAE algorithm [5] and very widely used by the engineering 
and architectural communities. The Iqbal model [8] offers extra-accuracy over more 
conventional models as reviewed by Gueymard [13]. 

Development of correlation models that predict the beam or sky radiation using 
other solar radiation measurements is possible. Decomposition models usually use 
information only on global radiation to predict the beam and sky components. 
These relationships are usually expressed in terms of the irradiations which are the 
time integrals (usually over 1 h) of the radiant flux or irradiance. Decomposition 
models developed to estimate direct and diffuse irradiance from global irradiance 
data were found in the literature [15-23]. 

Solar radiation models commonly used for systems using solar energy are 
reviewed in this article. Comparisons of predicted values with these models and 
measured data in Hong Kong are made. The units and symbols follow those sug¬ 
gested by Beckman et al. [24]. 


2. Parametric models 

2.1. Iqbal model C [8] 

The direct normal irradiance 7 n (W m -2 ) in model C described by Iqbal [8] is given by 

/„ = 0.9751 £o4; r r t 0 r g r w r a (1) 

where the factor 0.9751 is included because the spectral interval considered is 0.3-3 
pm; 7 SC is the solar constant which can be taken as 1367 W m -2 . E 0 (dimensionless) 
is the eccentricity correction-factor of the Earth’s orbit and is given by 

E a = 1.00011 + 0.034221cosr + 0.00128sinr + 0.000719cos2r + 0.000077sin2r 

(2) 

where the day angle T (radians) is given by 


where N is the day number of the year, ranging from 1 on 1 January to 365 on 31 
December. 
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T r , r 0 , T g , t w , r a (dimensionless) are the Rayleigh, ozone, gas, water and aerosols 
scattering-transmittances respectively. They are given by 

Tr = e -0.090 3m r(l+'»a-m‘ M ) (4) 


[o.1611E/ 3 (1 + 139.48C/ 3 r a3O35 -0.002715C/3(l + 0.044C/ 3 + 0.000317^) '1 


_ e -0.0127m»- 26 

r w - I - 2.4959C/i[(1 + 79.034t/ 1 ) 0(,828 + 6.3 8 5£/i]~ 

T = e -&,-(i + /„,-e° 8 K 9108 


where ra a (dimensionless) is the air mass at actual pressure and m T (dimensionless) is 
the air mass at standard pressure (1013.25 mbar). They are related by 


-'(isfe) 


(9) 


where p (mbar) is the local air-pressure. 

Ut, (cm) is the ozone’s relative optical-path length under the normal temperature 
and surface pressure (NTP) and is given by 


U 3 = lo7.m T 


( 10 ) 


where l oz (cm) is the vertical ozone-layer thickness. 

C7i (cm) is the pressure-corrected relative optical-path length of precipitable water, 
as given by 


U\ — wm r 


( 11 ) 


where w (cm) is the precipitable water-vapour thickness reduced to the standard 
pressure (1013.25 mbar) and at the temperature T of 273 K: w is calculated with the 
precipitable water-vapour thickness under the actual condition w' (cm) by 
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4o (dimensionless) is the aerosol optical thickness given by 

4o = 0.2758/ ao; ;.]A=0.38 nm + 0.354o,9.^=O.S urn (13) 

where X (pm) is the wavelength. 

The beam component 4, (W m -2 ) is given by 

7 b = cos0 z /„ (14) 

where 9 Z (degrees) is the zenith angle. 

The horizontal diffuse irradiance Id (W m -2 ) at ground level is a combination of 
three individual components corresponding to the Rayleigh scattering after the first 
pass through the atmosphere, D r (W m -2 ); the aerosols scattering after the first pass 
through the atmosphere, 7) a (W m -2 ); and the multiple-reflection processes between 
the ground and sky 7) m , (W m -2 ). 


Id = Dr + A, + An 


where 

A = 


0 .797 SC s i n a r„ r g r w r aa 0.5(1 — r r ) 
l-/ n a + ml 02 


(16) 


where a (degrees) is the solar altitude and is related to the zenith angle 9 Z by the 
following equation 

cos0 z = si no; (17) 


r aa (dimensionless) is the transmittance of direct radiation due to aerosol absorp- 
tance and is given by 


*aa = 1 - (1 - «o)(l - m a + 777 ‘ 06 )(1 - T a ) (18) 

where co 0 (dimensionless) is the single-scattering albedo fraction of incident energy 
scattered to total attenuation by aerosols and is taken to be 0.9. 

n 0.797 SC s i n a r 0 r g r w r aa T c (1 - r as ) /1fu 

a _ 1 - 7u a + mi- 02 ( } 

F c (dimensionless) is the fraction of forward scattering to total scattering and is 
taken to be 0.84; r as (dimensionless) is the fraction of the incident energy transmitted 
after the scattering effects of the aerosols and is given by 


( 20 ) 
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D m 


(/„sina + £>i + £>a)p g Ai 

1 - PgPa 


( 21 ) 


where p g (dimensionless) is the ground albedo; p a (dimensionless) is the albedo of 
the cloudless sky and can be computed from 


p a = 0.0685 + (1 - F c ) (1 - r as ) (22) 

where (1 — F c ) is the back-scatterance. Consequently, the second term on the right- 
hand side of this equation represents the albedo of cloudless skies due to the pre¬ 
sence of aerosols, whereas the first term represents the albedo of clean air. 

The global (direct plus diffuse) irradiance 7 t (W m -2 ) on a horizontal surface can 
be written as 

7, = h + h = (/nCosPz + 7) r + £> a ) ) (23) 

2.2. ASHRAE model [5] 

A simpler procedure for solar-radiation evaluation is adopted in ASHRAE [5] and 
widely used by the engineering and architectural communities. The direct normal 
irradiance 7 n (W m -2 ) is given by 

7 n = (C n )Ae-i^ (24) 

where A (W m -2 ) is the apparent extraterrestrial irradiance, as given in Table 1, and 
which takes into account the variations in the Sun-Earth distance. Therefore 


Table 1 

Values of A, B and C for the calculation of solar irradiance according to the ASHRAE handbook: 
applications [5] 


Date 

A (W m- 2 ) 

B (W in A 

C (W m 2 ) 

21 Jan 

1230 

0.142 

0.058 

21 Feb 

1215 

0.144 

0.060 

21 Mar 

1186 

0.156 

0.071 

21 Apr 

1136 

0.180 

0.097 

21 May 

1104 

0.196 

0.121 

21 Jun 

1088 

0.205 

0.134 

21 Jul 

1085 

0.207 

0.136 

21 Aug 

1107 

0.201 

0.122 

21 Sep 

1152 

0.177 

0.092 

21 Oct 

1193 

0.160 

0.073 

21 Nov 

1221 

0.149 

0.063 

21 Dec 

1234 

0.142 

0.057 
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A — 7 sc 7i 0 (constant) (25) 

The value of the multiplying constant is close to 0.9 [8]. The variable B (dimen¬ 
sionless) in Table 1 represents an overall broadband value of the atmospheric 
attenuation coefficient for the basic atmosphere of Threlkeld and Jordan [5]. C n 
(dimensionless) is the clearness number and the map of C n values for the USA is 
provided in the ASHRAE handbook: applications [5]. C n is the ratio of the direct 
normal irradiance calculated with the local mean clear-day water-vapour to the 
direct normal irradiance calculated with water vapour according to the basic atmo¬ 
sphere. 

Eq. (24) was developed for sea-level conditions. It can be adopted for other 
atmospheric pressures by 

/„ = (C n )Ae ”*Go) = C n Ae _jB fe) sec(?z (26) 

where p (mbar) is the actual local-air pressure and p Q is the standard pressure 
(1013.25 mbar). 

In the above equation, the term (p/p 0 ) sec 0 7 approximates to the air mass, with the 
assumptions that the curvature of the Earth and the refraction of air are negligible. 

The total solar irradiance 7 t (W m -2 ) of a terrestrial surface of any orientation and 
tilt with an incident angle 9 (degrees) is the sum of the direct component 7 n cos 9 plus 
the diffuse component 7 d coming from the sky, plus whatever amount of reflected 
short-wave radiation 7 r (W m -2 ) that may reach the surface from the Earth or from 
adjacent surfaces, i.e., 


I t = 7 n cos6> + 7 d + 7 r (27) 

The diffuse component is defined through a variable C (dimensionless) given in 
Table 1. The variable C is the ratio of the diffuse radiation falling on a horizontal 
surface under a cloudless sky 7 d to the direct normal irradiation at the Earth’s sur¬ 
face on a clear day 7 n . The diffuse radiation 7 d is given by 


7 d = CI n F 5 S 


(28) 


where for a horizontal surface, 



(29) 


F ss (dimensionless), the angle factor between the surface and the sky, is given by 


F ss 


1 + cosf 
2 


( 30 ) 


where f (degrees) is the tilt angle of a surface measured from the horizontal. 
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The reflected radiation 7 r (W m 2 ) from the foreground is given by 

K — iv.0= 0 K g F sg (31) 

where /c g (dimensionless) is the reflectance of the foreground; and / t; e=o is the total 
horizontal irradiation (9 — 0°). F sg , the angle factor between the surface and the 
earth, is given by 


_ _ i - cost: 

<sg = ^— 


(32) 


3. Comparison of the Iqbal model C [8] and ASHRAE model [5] 

The direct normal irradiance I n was calculated with the Iqbal model C [8] for an 
ozone thickness of 0.35 cm (NTP) at the 21st day of each month for various values 
of the zenith angle 0 Z for Hong Kong (22.3°N latitude, 114.3°E longitude). The 
beam, diffuse and global irradiances 7b, id and 7 t respectively, were computed using 
both models and are shown in Fig. l(a)-(c). The predicted values for January and 
June are shown in Fig. 2(a)-(c) for comparison. 

At 6 Z — 0°, a maximum difference for the beam irradiance of not more than 7% is 
obtained between the two models. The beam irradiance deviation of less than 8% was 
found for September-April, and less than 12% for May-August when 0 7 < 60°. The 
ASHRAE model would be inaccurate for the following reasons [2,5]. First, the solar 
constant used (1322 W m -2 ) and the extraterrestrial spectral irradiance employed is 
outdated. Secondly, the attenuation coefficients are empirical functions and are based on 
site-specific data (Mount Wilson and Washington), and old instruments were used to 
take the data. Thirdly, when (In /„) versus (p/p 0 ) sec 0 7 was plotted from Eq. (26), the 
result was a straight line. The ordinate gives the apparent solar constant A, and the 
slope of the straight line is the overall broadband atmospheric attenuation-coefficient 
B. This linear relationship would be the approximation for broadband irradiance. 

A large deviation (up to 41%) was found between the predicted values of diffuse 
irradiance with the two models as shown in the figure. For the ASHRAE model, the 
ground albedo is not included, so contributions by multiple reflections are ignored. 
The aerosol-generated diffuse irradiance is not considered in the ASHRAE model. 
This may introduce a large error as shown in Figs. 1(b) and 2(b). However, as dif¬ 
fuse radiation is only a small fraction of the direct radiation, the error is not serious 
in many applications. The deviation of the global irradiance 7 t from the predicted 7 t 
at 6 Z at 0° for all months was only 5%. The predicted 7 t shown in Figs. 1(c) and 2(c). 


4. Correlation of average daily global radiation with hours of sunshine 

Theoretical determinations of the direct, diffuse and directional intensities of dif¬ 
fuse irradiance would require data on the type and optical properties of clouds, 
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cloud amount, thickness, position and the number of layers. These data are rarely 
collected on a routine basis. However, sunshine hours and total cloud-cover data 
(i.e. the fraction of sky covered by clouds) are widely and easily available. Correla¬ 
tions, therefore, were developed to estimate insolation and sunshine hours. 

The monthly average daily global radiation on a horizontal surface can be esti¬ 
mated through the number of bright sunshine hours. Angstrom [25] developed the 
first model suggesting that the ratio of the average daily global radiation G t (MJ 
m -2 ) and cloudless radiation G c (MJ m -2 ) is related to the monthly mean daily 
fraction of possible sunshine (sunshine fraction) Sf (dimensionless) by 


Gt 

Gc 


— a\ + (1 — ui)*Sf — a\+ biSf 


(33) 


with the constant a\ being obtained as 0.25 for Stockholm, Sweden. 

The sunshine fraction Sf is obtained from: 

St = §- o ( 34 ) 

where S (h) is the monthly average number of instrument-recorded bright daily 
sunshine hours, and S Q (h) is the average day length. In the Angstrom model, 
a\ + b\ — \, because on clear days, S { is supposed to equal unity. However, because 
of problems inherent in the sunshine recorders, measurements of S f would never 
equal unity. This leads to an underestimation of the insolation. 

A modified form of the above equation, known as the Angstrom-Prescott for¬ 
mula, was reported by Prescott in 1940 [26]: 


-Q - — <2i + b\S{ 


(35) 


where a\ = 0.22 and b\ = 0.54. 

In Eqs. (33) and (35), a\ and b\ correspond respectively to the relative diffuse 
radiation on an overcast day, whereas ( a\ +b\) corresponds to the relative cloud¬ 
less-sky global irradiation. These equations assume that G t corresponds to an idea¬ 
lised day over the month and can be obtained as the weighted sum of the global 
irradiations accumulated during two limiting days with extreme sky conditions, i.e. 
fully overcast and perfectly cloudless. 

There may be problems in calculating G c accurately. This equation was modified 
to be based on the extraterrestrial irradiation by Prescott in 1940 as reviewed by 
Gueymard et al. [26]: 


G t , S 

G=“' +b 'Y„ 


( 36 ) 
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In equation (36), G t and G a (MJ m -2 ) are the monthly mean daily global and 
extraterrestrial radiations on a horizontal surface, ci\ and b\ are constants for dif¬ 
ferent locations. 

The hourly extraterrestrial radiation on a horizontal surface G 0 (MJ m -2 ) can be 
determined by [8] 

G a — — / sc £ 0 cos</>cos<5^sin&j s — (]^)joj s cosoj s j (37) 

where 7 SC is the solar constant (=1367x3.6 kJ m -2 h -1 ), and co s (degrees) is the 
sunset-hour angle for a horizontal surface. 

Rietveld [27] examined several published values of a x and b\ and noted that a t is 
related linearly and b\ hyperbolically to the appropriate yearly average value of 5 f , 


denoted as (J^), such that 


ai =0.10 + 0.24 

(38) 

© 

© 

+ 

ro 

© 

II 

(39) 

If the Rietveld model [27] can be simplified [26] to a constant-coefficient 

Angstrom-Prescott equation by substituting Eqs. (38) and (39) into (36): 

^ = 0.18 + 0.62^ 

(40) 

This equation was proposed for all places in the world and 
superior results for cloudy conditions when < 0.4. 

Glover and McCulloch included the latitude <]) (degrees) effect 
following correlation, as reviewed by Iqbal [8]: 

yields particularly 

and presented the 

^ = O.29cos0 + 0.52-^- for <jx 60° 

G 0 o 0 

(41) 

Gopinathan [28] suggested the regression coefficients a x and b\ in terms of the 
latitude, elevation and percentage of possible sunshine for any location around the 
world. The correlations are: 

a x = -0.309 + O.539cos0 - 0.0693 h + 0.29 

(42) 

bx = 1.527 - l.O27cos0 + 0.09266 - 0.359^ 

*^0 

(43) 
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where 0 (degrees) is the latitude, and h (km) is the elevation of the location above sea 
level. 

The coefficients ci\ and b\ are site dependent [20,29-32]. A pertinent summary is 
shown in Table 2. The correlations of the average daily global radiation ratio and 
the sunshine duration ratio in Hong Kong (22.3°N latitude) are shown in Fig. 3 for 
comparison. The monthly average daily global radiation ratio registered in 1961 
1990 in Hong Kong [33] is shown. It is very close to the predictions of Bahel et al. 
[29], Alnaser [30] and Louche et al. [20], 

These coefficients would be affected by the optical properties of the cloud cover, 
ground reflectivity, and average air mass. Hay [34] studied these factors with data 
collected in western Canada and proposed a site-independent correlation 


Gt 

Go 


0.1572 + 0.556 



(44) 


where p g (dimensionless) is the ground albedo, p a is the cloudless-sky albedo taken 
as 0.25 and p c is the cloud albedo taken as 0.6. »S om (h) is the modified day-length for 
a solar zenith angle 9 Z greater than 85° and is given by 


1 _i /cos85 — sin</)sin<5 c \ 
7.5 \ cos0cos<5 c / 


(45) 


where <5 C (degrees) is the characteristic’s declination, i.e. the declination at which the 
extraterrestrial irradiation is identical to its monthly average value. 

Sen [32] employed fuzzy modelling to represent the relations between solar irra¬ 
diation and sunshine duration by a set of fuzzy rules. A fuzzy logic algorithm has 
been devised for estimating the solar irradiation from sunshine duration measure¬ 
ments. The classical Angstrom correlations are replaced by a set of fuzzy-rule bases 


Table 2 

Some examples of regression coefficients a\ and b t used in different models 


Reference 

ai 

bi 

Data source 

[27] 

0.1+0.24 ( 1 ) 

0.38 + 0.08(f) 

World-wide (42 places) (6° < 0 < 69°) 

Simplified [27] 

0.18 ' o/ 

0.62 ' o/ 

World-wide (42 places) (6° < 0 < 69°) 

[8] 

O.29cos0 

0.52 

0 < 60° 

[29] 

0.175 

0.552 

Dhahran, Saudi Arabia 
(17° < 0 < 27°) 

[28] 

-0.309+ 0.539cos 0- 
0.0693h + 0.29 f 

1.527—l.O27cos0 + 
0.0926h—0.3591 

World-wide (5° <0<54°) 

[30] 

0.2843 

0.4509 

Bahrain (0 = 27°) 

[ 20 ] 

0.206 

0.546 

Southern France (0 = 42°) 

[31] 

0.2006 

0.5313 

Northern India (0 = 27°) 
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Sunshine duration ratio S/S 0 

Fig. 3. Correlation of average daily global radiation with sunshine duration for Hong Kong (t/> = 23°). 


and applied for three sites with monthly averages of daily irradiances in the western 
part of Turkey. 

Mohandes et al. [35] studied the monthly mean daily values of solar radiation 
falling on horizontal surfaces using the radial basis functions technique. The pre¬ 
dicted results by a multilayer perceptrons network were compared with the simpli¬ 
fied Rietveld model [27]. The predicted results by the neural network were compared 
with those measured data for ten locations and good agreement was found. 


5. Estimation of monthly-average hourly global radiation 


As reviewed by Iqbal [8], Whiller investigated the ratios of hourly global radiation 
7 t (MJ m -2 ) to daily global radiation G t (MJ m -2 ) in 1956 and assumed that 


I t _ I 0 
G~ t ~ Go 


( 46 ) 
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where I 0 is the extraterrestrial hourly radiation (MJ m 2 ), which can be deter¬ 
mined by 

I 0 = / Sc 7io(sin<5sin0 + O.9972cos<$cos0coso>i) (47) 

where I sc is the solar constant, E 0 is the eccentricity correction factor, S, <p and 0\ are 
the declination, latitude and hour angle respectively at the middle of an hour in 
degrees. 

Whiller [8] correlated the hourly and daily beam radiations 7b (MJ m -2 ) and Gb 
(MJ m -2 ) respectively, as 

h_ _ Io_ _ S Xh k h (w)i x E 0 cose z &co , 

Gb G 0 J 1 day /c b (&j)/ sc G o cos0 z doj 


where w (degrees) is the hour angle, i.e., zero for solar noon and positive for the 
morning. 

Assuming that the atmospheric transmittance for beam radiation kb(a >) to be 
constant throughout the day, this equation can be rewritten as 


Ib_ _ Io_ 

Gb G 0 


24 . (ix\ 

„ —sin — Icoscui - cosojc 
71 it V24/ _ 

24 sinews - (y^joJsCOSOJs 


(49) 


where w s is the sunset-hour angle (in degrees) for a horizontal surface. 

Liu and Jordan [15] extended the day-length range of Whiller’s data and plotted 
the ratio >\ = ^ against the sunset-hour angle &j s . The mathematical expression was 
presented by Collares-Pereira and Rabl [36] 


n 


I t 

G~t 


(50) 


where 

a 2 = 0.409 + 0.5016sin(w s - 60°) 
b 2 = 0.6609 - 0.4767sin(w s - 60°) 


With Eqs. (37) and (47) for I 0 and G 0 , Eq. (50) can be expressed as 


n 


71 

24 



(a 2 + bjcosojf) 


(51) 

(52) 


(53) 
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A graphical presentation is shown in Fig. 4 and the graph is applicable everywhere 
[e.g. 8,15,29]. 


6. Estimation of the hourly diffuse radiation on a horizontal surface using 
decomposition models 

Values of global and diffuse radiations for individual hours are essential for 
research and engineering applications. Hourly global radiations on horizontal sur¬ 
faces are available for many stations, but relatively few stations measure the hourly 
diffuse radiation. Decomposition models have, therefore, been developed to predict 
the diffuse radiation using the measured global data. 

The models are based on the correlations between the clearness index k t (dimen¬ 
sionless) and the diffuse fraction k<\ (dimensionless), diffuse coefficient ku (dimen¬ 
sionless) or the direct transmittance kb (dimensionless) where 



Sunset hour angle w s / degrees 


8 9 10 11 12 13 14 15 16 

Hours from sunrise to sunset 


Fig. 4. Distribution of the monthly average hourly global radiation. 
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k t 


h 

4 


(54) 


k d 


Id 

/, 


ku — 


h 

Io 


(55) 

(56) 


k h 


4 

4 


(57) 


It, 4, 4 and 4 being the global, direct, diffuse and extraterrestrial irradiances, 
respectively, on a horizontal surface (all in MJ m -2 ). 

6.1. Liu and Jordan model [15] 

The relationships permitting the determination, for a horizontal surface, of the 
instantaneous intensity of diffuse radiation on clear days, the long-term average 
hourly and daily sums of diffuse radiation, and the daily sums of diffuse radiation 
for various categories of days of differing degrees of cloudiness, with data from 98 
localities in the USA and Canada (19° to 55°N latitude), were studied by Liu and 
Jordan [15], 

The transmission coefficient for total radiation on a horizontal surface is given by 
the intensity of total radiation (i.e. direct 4 plus diffuse 4) incident upon a hor¬ 
izontal surface 7 t divided by the intensity of solar radiation incident upon a hor¬ 
izontal surface outside the atmosphere of the Earth I Q . The correlation between the 
intensities of direct and total radiations on clear days is given by 

ku — 0.271 — 0.29394, (58) 


Since 


k t = (7b + 7n)/4 = 4> + k n 


(59) 


then 


ku = 0.384 - 0.4164 


(60) 


6.2. Orgill and Hollands model [16] 

Following the work by Liu and Jordan [15], the diffuse fraction was estimated by 
Orgill and Hollands [16] using the clearness index 4 as the only variable. The model 
was based on the global and diffuse irradiance values registered in Toronto (Canada, 
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42.8°N) in the years 1967-1971. The correlation between hourly diffuse fraction k d 
and clearness index k t is given by 


k d = 1 - 0.249k t ; k, < 0.35 

(61) 

k d = 1.577 - 1.84k,; 0.35 <k,< 0.75 

(62) 

k d = 0.177; k,> 0.75 

(63) 


Once the hourly diffuse fraction is obtained from the k d —k t correlations, the 
direct irradiance 7b is obtained by [23] 


St(l~*d) 

sinu 


(64) 


where a is the solar elevation angle. 

6.3. Erbs et al. model [17] 

Erbs et al. [17] studied the same kind of correlations with data from 5 stations in 
the USA with latitudes between 31 and 42°. The data were of short duration, ran¬ 
ging from 1 to 4 years. In each station, hourly values of normal direct irradiance and 
global irradiance on a horizontal surface were registered. Diffuse irradiance was 
obtained as the difference of these quantities. The diffuse fraction k d is calculated by 

k d = 1 - 0.09k, for k,< 0.22 (65) 

k d = 0.9511 -0.1604k,+ 4.388^+ 16.638^+ 12.336^ for 0.22<k,<0.8 

(66) 

k d = 0.165 for k t >0.8 (67) 


6.4. Spencer model [21] 

Spencer [21] studied the latitude dependence on the mean daily diffuse radiation 
with data from 5 stations in Australia (20^45°S latitude). The following correlation 
is proposed 


k d = fl 3 - hki for 0.35 <k t < 0.75 ( 68 ) 

It is presumed that k d has constant values beyond the above range of k t . The 
coefficients a 3 and 63 are latitude dependent. 


« 3 = 0.94 + 0.0118101 


(69) 
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*3 = 1.185 + 0.0135101 (70) 

where 0 (degrees) is the latitude. There would be an increasing proportion of diffuse 
radiation at higher latitudes due to the increase in average air mass. 

6.5. Reindletal. [19] 

Reindl et al. [19] estimated the diffuse fraction k d using two different models 
developed with measurements of global and diffuse irradiance on a horizontal sur¬ 
face registered at 5 locations in the USA and Europe (28-60°N latitude). The first 
model (denoted as Reindl-1) estimates the diffuse fraction using the clearness index 
as input data and the diffuse fraction * d is given by 

* d = 1.02 - 0.248frt for *, < 0.3 (71) 

k A = 1.45 - 1.67*, for 0.3 < *, < 0.78 (72) 

* d = 0.147 for *,>0.78 (73) 


The second correlation (denoted as the Reindl-2 model) estimates the diffuse 
fraction in terms of the clearness index *, and the solar elevation a. The diffuse 
fraction * d is given by 

* d = 1.02 — 0.254*, + 0.0123sina for *,<0.3 (74) 

* d = 1.4 - 1.749*, + 0.177sina for 0.3 < *, < 0.78 (75) 

* d = 0.486*, - 0.182sina for *, ^0.78 (76) 


6.6. Lam and Li [22] 

Lam and Li [22] studied the correlation between global solar radiation and its 
direct and diffuse components for Hong Kong (22.3°N latitude) with the measured 
data in 1991-1994. A hybrid correlation model based on hourly measured data for 
the prediction of hourly direct and diffuse components from the global radiation for 
Hong Kong is given by 

* d = 0.977 for *,<0.15 (77) 

* d = 1.237 - 1.361*, for 0.15 < *,<0.7 (78) 


* d = 0.273 for *, > 0.7 


(79) 
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6.7. Skartveit and Olseth model [18] 

Skartveit and Olseth [18] showed that the diffuse fraction depends also on other 
parameters such as solar elevation, temperature and relative humidity. Similar 


arguments were found in the literature [19,37-39]. 

Skartveit and Olseth [18] estimated the direct irradiance 7 b from the global irrad- 
iance G t and from the solar elevation angle a for Bergen, Norway (60.4°N latitude) 
with the following equation 

r G t (l - f) 

•*b — 

sina 

(80) 

where x] is a function of k t and the solar elevation a (degrees). The model was vali¬ 
dated with data collected in Aas, Norway (59.7°N latitude), Vancouver, Canada 
(49.3°N latitude) and 10 other stations worldwide. Details of this function are given 
below: 

If A t < Ci, 


(81) 

where c\ — 0.2. 

If ci <k t < 1.09 c 2 , 


+ = 1 - (1 - di)[d 2 cl /2 + (1 - d 2 )c^ 

(82) 

where 


c 2 = 0.87 - 0.56e _0 06 “ 

(83) 


(84) 

C4 — k t — Cl 

(85) 

di =0.15 + 0.43e _006 “ 

(86) 

d 2 = 0.27 

(87) 

d 2 = c 2 - ci 

(88) 

If At > 1.09 c 2 



yjr= 1 - 1.09c 2 V^ 


(89) 
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where 


l-(l-d l )[d2C / 3 - + (l-d 2 y 3 2 J 

(90) 

K i+ 4(!-°- 5 )]) 

(91) 

1.09c 2 - ci 

(92) 


The constants have to be adjusted for conditions deviating from the validation 
domain. 

6.8. Maxwell model [23] 

A quasi-physical model for converting hourly global horizontal to direct normal 
insolation proposed by Maxwell in 1987, was reviewed by Bathes et al. (23). The 
model combines a clear physical model with experimental fits for other conditions. 
The direct irradiance 7 b is given by: 

/b = /o{V'-(^4 + ^c m “‘ 4 )} (93) 

where I a is the extraterrestrial irradiance, is a function of the air mass m a 
(dimensionless) and is given by 

= 0.866 - 0.122m a + 0.0121m 3 - 0.000653m 3 + 0.000014m|| (94) 

and d 4 , d 5 and d 6 are functions of the clearness index k t as given below: 

IfA t <0.6, 


d 4 = 0.512 - 1.56 k t + 2.286 *£ - 2.222 k 3 
d 5 = 0.37 + 0.962 k t 

d 6 = - 0.28 + 0.923£ t —2.048 k\ (95) 

If At > 0.6, 

<k = -5.743 + 21.77A, - 27.49^ + 11.56^ 
d s = 41.4 - 118.5At + 66.05 k\ + 31.9A? 

d 6 =t -47.01 + 184.2A, - 222^ + 73.81A t 3 (96) 


6.9. Louche et al. model [20] 

Louche et al. [20] used the clearness index k t to estimate the transmittance of beam 
radiation k b . The correlation is given by 
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* b = —10.627 k 5 t + 15.307 k\ - 5.205 k\ + 0.994*? - 0.059*, + 0.002 (97) 

The correlation includes global and direct irradiance data for Ajaccio (Corsica, 
France, 44.9°N latitude) between October 1983 and June 1985. 

Other researchers [40,41] estimated the direct irradiance by means of k\>—k t cor¬ 
relations. They found that the solar elevation is an important variable in this type of 
correlation. The direct irradiance is then estimated with the following definition of 
direct transmittance [23] 

h = k h I 0 (98) 


6.10. Vignola and McDaniels model [40] 

Vignola and McDaniels [40] studied the daily, 10-day and monthly average beam- 
global correlations for 7 sites in Oregon and Idaho, USA (38^46°N latitude). The 
beam-global correlations vary with time of year in a manner similar to the seasonal 
variations exhibited by diffuse-global correlations. The direct transmittance *b for 
daily correlation is given by 

k h = 0.022 - 0.28*, + 0.828*2 + 0.765*? for *, 2*0.175 (99) 

* b = 0.016*, for *, < 0.175 (100) 

and, with seasonal terms, 


* d = 0.013 - 0.175*, + 0.52*2 + , mk i 

+ (0.038*, - 0.13*2)sin[27r(A- 20)/365] (101) 

for *,>0.175 

* b = 0.125*? for *, < 0.175 (102) 

where N is the year-day under consideration. 

6.11. Al-Riahi et al. model [42] 

Al-Riahi et al. [42] developed a clear-day model for the beam transmittance of the 
atmosphere based on 5-year daily global radiation data measured at the Solar 
Energy Research Center (Baghdad, Iraq, 33.2°N latitude). The model offers the 
prediction of clear-day hourly values of direct-normal and global solar radiations 
for any day of the year at a given location with no required meteorological inputs. 
The global irradiance /, is given by 
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/, = / b (sina + f) (103) 

where ^ is a fraction of the diffuse irradiance divided by the direct irradiance and 
would be a constant throughout a particular clear-day. The fraction xpata given day 
N (dimensionless) of a year is given by 

x!f(N) = 0.0936 + 0.041sin p iV ~ i 1 6 0 7 4 ' 5)7r j + 0.004773sin|^^±|l^J (i 0 4) 


4 is given by 

4 = 4e -ft " a = 4 e- jB (as)° 678 (105) 


where m a is the air mass, and B is the overall atmospheric coefficient of extinction 
which varies over the year. The exponent 0.678 compensates for the actual curvature 
of the atmosphere. 

B(N) is given by 

B(N) = 0.3917397 - 5.596 x 10_2sin (^^ + 5.293 x 10_3cos (^ 7V ) 

+ 1.3594 x l°" 2 sin^^ +4.0383 x 10_3cos (^^ 

(106) 


7. Comparison of decomposition models 

The k a —k t correlations of Liu and Jordan [15], Orgill and Hollands [16], Erbs et 
al. [17], Reindl-1 [19], Spencer [21] and Lam and Li [22] are plotted in Fig. 5. The 
correlations of Orgill and Hollands [16], Erbs et al. [17], Reindl-1 [19] and Lam and 
Li [22] are almost identical. The predicted diffuse fraction for Hong Kong by Lam 
and Li’s correlations is higher than by other correlations for clear-sky conditions 
(k t > 0.7). The Spencer [21] correlation plotted for (p yields consistently lower 
results. 

The kd — k t correlations of Reindl-2 [19] and Skartveit and Olseth [18] were plot¬ 
ted in Fig. 6 for solar elevations of 10, 40 and 90°. The results agree well with each 
other for k t ^0.3 and agree reasonably well for 0.3 < k t < 0.8. A larger deviation 
was found for k t > 0.8 and for a lower solar angle a — 10°. The predicted diffuse 
fraction for Hong Kong [22] was plotted in the figure for comparison. The predicted 
results agree with those of the other two models. 

The kb— k t correlations proposed by Louche et al. [20], Vignola and McDaniels 
[40] and Maxwell [23] with air masses of 1, 2 and 4 are plotted in Fig. 7. The results 
of Louche et al. [20] and Vignola and McDaniels [40] are almost identical. The 
Maxwell modelled results agree well with the predictions of the other two models. The 
kb — k t correlation for Hong Kong would be derived from the kd—k t correlation 
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Clearness index k t 

Fig. 5. Modelled hourly diffuse fraction k d as a function of clearness index k t . 


[Eqs. (77) to (79)] proposed by Lam and Li [22], With Eqs. (54), (55) and (57), k b can 
be expressed as 


k b - (1 - k A )k t 


(107) 


The derived k b —k t correlation for Hong Kong agrees well with k b predicted by the 
others for k t < 0.7 and falls below the predictions for k t > 0.7 as shown in Eig. 7. 


8. Comparison of decomposition models and parametric model 

Batlles et al. [23] estimated the hourly values of direct irradiance using the 
decomposition models [16-20]. The results of these models were compared with 
those provided by an atmospheric transmittance parametric model [8] and with the 
measured direct irradiances for 6 stations (36.8 43°N latitude, 1.7-5.4°W longitude) 
reported by Batlles et al. [43]. 

The aerosol-scattering transmittance adopted by Batlles et al. [23] is calculated 
from the visibility. 










L.T. Wong, W.K. Chow / Applied Energy 69 (2001) 191-224 


Ta = (0.97 - 1.265F^°' 66 ) m “ 9 (108) 

where visibility F; s (km) [2] is given by 

Vi s — 147.994- 1740.523^-(£> 2 -0.17^ +0.011758) 0 ' 5 ] (109) 


and 


yjr = 0.55" a (110) 

where fix and 7*2 are the Angstrom turbidity-parameters: P 2 is taken to be 1.3 and fix 
is calculated using the Linke turbidity-factor 

71 = ( 39.5e-- :> '+47.4 + t ’' 1 ) +(16 + ° 22 ” /)ft (111) 

where a is the solar elevation in degrees and w' is the precipitable water-content in 
cm. The Linke turbidity-factor T L is defined as the number of Rayleigh atmospheres 
(an atmosphere clear of aerosols and without water vapour) required to produce a 
determined attenuation of direct radiation and is expressed by 

7T = (/R 0 m a r'h/ 0 (112) 

lb 


where / Ro (dimensionless) is the Rayleigh optical-thickness. In the study by Bathes et 
al. [23], hourly and monthly average hourly values of were used. 

The decompositions are good for high values of the clearness index (i.e. cloudless 
skies and high solar elevations). In such conditions, the model proposed by Louche 
et al. [20] estimated the direct irradiance with a 10% root-mean-square (RMS) error 
and the estimated direct irradiance by the Iqbal model [8] using hourly values of fi 
has a 4% RMS error. Both models have no significant mean bias error. Bathes et al. 
[23] concluded that the parametric model is the best when precise information of the 
turbidity coefficient is available. However, if there is no turbidity information 
available, decomposition models are a good choice. 

The clearness index k t and diffuse fraction k& for Hong Kong (22.3°N latitude) were 
estimated with the parametric model C [8] as this model offers accurate prediction 
over a wide range of conditions [13]. The aerosols, Rayleigh, ozone, gas and water 
scattering transmittances r a , r r , r 0 , r g and r w respectively, were determined by Eqs. 
(108) and (4)-(7). The Angstrom turbidity parameter fij is taken to be 1.3 and fi\ is 
taken to be 0, 0.1, 0.2 and 0.4 for the clean, clear, turbid and very turbid atmo¬ 
spheres respectively [8]. The corresponding visibility Vi s in Eq. (109) is 337, 28, 11 
and 5 km, respectively. The measured vertical ozone-layer thickness l oz for the 
Northern Hemisphere [44] was used to determine the ozone relative optical-path 
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length U 3 in Eq. (10). The air mass at standard pressure, m r , is approximated by 
Kasten’s formula [8] 

m r = [cosd z + 0.15(93.885 - dz) -1 ' 253 ]” 1 (113) 

The precipitable water-vapour thickness W is approximated by Won’s equation [8] 

yj = 0. le 2.2572+0.05454r dcw (j 14 ) 

where T dcw is the dew-point temperature (°C). The monthly average daily tempera¬ 
ture T and monthly average daily dew-point temperature T dcv/ recorded between 
1961 and 1990 in Hong Kong [33] were adopted in the calculations. The local air 
pressure p is taken to be 1013.25 mbar. The ground albedo p g is taken to be 0.09, 
0.35 and 0.74 for black concrete, uncoloured concrete and white glazed surfaces 
respectively (as Hong Kong is a developed city of many high-rise concrete buildings 
with curtain walls). 

The extraterrestrial radiation 7 0 was determined by [8] 

7 0 = 7 sc £ 0 cos6» z (115) 

The clearness index k t and the diffuse fraction k d were determined with the mod¬ 
elled global solar irradiance 7 t and total diffuse irradiance 7 d . The average values of 
the predicted diffuse fraction by the Iqbal model C [8] for clean, clear, turbid and 
very turbid atmospheres at ground albedo of 0.09, 0.35 and 0.74 are plotted in Fig. 8. 
The predicted k d by Lam and Li’s correlations [22] for Hong Kong is shown for 
comparison. Lam and Li’s results agree with the modeled results for F; s of about 28 
km for k t between 0.4 and 0.7 and fall in the range of Vi S between 11 and 28 km for 
k t between 0.2 and 0.4. 


9. Other studies of solar radiation models 

Literature results [37,38] showed that the diffuse fraction depends also on other 
variables like atmospheric turbidity, surface albedo and atmospheric precipitable 
water. 

Chendo and Maduekwe [39] studied the influence of four climatic parameters on the 
hourly diffuse fraction in Lagos, Nigeria. It was found that the diffuse fraction of the 
global solar-radiation depended on the clearness index k t , the ambient temperature, 
the relative humidity and the solar altitude. The standard error of the Liu and Jordan- 
type [16,17,19] was reduced by 12.8% when the solar elevation, ambient temperature 
and relative humidity were included as predictor variables for the 2-year data set. 
Tests with existing correlations showed that the Orgill and Hollands model [16] 
performed better than the other two models [17,19]. 

Elagib et al. [45] proposed a model to compute the monthly average daily global 
solar-radiation G t (MJ m -2 ) for Bahrain (22°N latitude) from commonly measured 
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Clearness index k t 

Fig. 8. Modelled hourly diffuse fraction k a as a function of clearness index k t for Hong Kong. 


meteorological parameters. The parameters are monthly average daily relative- 
humidity RH (%), the monthly average daily temperature-range AT (°C), the 
monthly average daily temperature T (°C), the monthly average daily hours of sun¬ 
shine S (h) and the monthly average daily extraterrestrial solar-radiation G c (MJ 
m“ 2 ). 

For January-June, G t can be computed by one of the following correlation equations: 


G t = 70.8385 - 0.7886 RH (116) 

G t = 55.8528 - 0.6440 (RH - S) (117) 

G t = 49.2637 — 0.6751(7?77— AT — S) (118) 

For July-December, 

G t = 4.5456e 00453r (119) 
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G t — —1.5358 — 0.6078G o (120) 

G t = 34.5577 - 0.4642 (RH - T) (121) 

G t = 18.9469 - 0.2663 (RH - T-G 0 ) (122) 


The 3-variable models yielded the smallest errors yet the 1-variable models can 
give close estimations. 

Tovar et al. [46] studied the probability density distributions of 1-min values of 
global irradiance, conditioned to the optical air-mass, considering them to be an 
approximation to the instantaneous distributions. It was found that the bimodality that 
characterizes these distributions increases with optical air-mass. The function based on 
Boltzmann’s statistics can be used for the generation of synthetic radiation data. 
Expressing the distribution as a sum of two functions provides an appropriate mod¬ 
elling of the bimodality feature that can be associated with the existing two levels of 
irradiance corresponding to two extreme atmospheric situations — the cloudless and 
cloudy conditions. 

In addition to the shortwave (0.3-2.0 pm) radiation from the Sun, the Earth 
receives longwave radiation (4—100 pm) from the atmosphere. The air conditions at 
ground level largely determine the magnitude of the incoming radiation [5] and the 
downward radiation from the atmosphere q iat could be expressed as 

7rat = £-dt<7 T^t (123) 


where e at (dimensionless) is the atmosphere emittance, a is the Stefan-Boltzmann 
constant (5.67xl0 -8 W m -2 K -4 ) and T. dt (K) is the ground-level temperature. 

The apparent sky-temperature T s ky (K) is defined as the temperature at which the 
sky (as a blackbody) emits radiation at the rate actually emitted by the atmosphere 
at ground-level temperature with its actual emittance e at . Then, 

<77^ = e at <r7< (124) 

or 

2"sky = GuTtt (125) 

The atmosphere emittance is a complex function of air temperature and moisture 
content [47]. A simple relationship [5], which ignores the vapour pressure of the 
atmosphere could be used to estimate the apparent sky-temperature 

T sky = 0.05527^ 5 


(126) 
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10. Conclusion 

Solar radiation models are desirable for designing solar-energy systems and good 
evaluations of thermal environments in buildings. Two parametric models (Iqbal 
model C [8] and ASHRAE model [5]) were used to predict the beam, diffuse and 
global irradiances under clear-sky conditions for Hong Kong. At zero zenith angle, 
a maximum difference of beam irradiance of not more than 7% is obtained between 
the 2 models. The beam irradiance deviation of less than 8% was found for Sep- 
tember-April, and less than 12% for May-August when the zenith angle is less than 
60°. The deviation of the global irradiance predicted between the 2 models was only 
5%. A large deviation, up to 41%, was found between the predicted values of the 
diffuse irradiance with the 2 models. The ASHRAE model would be less accurate 
for diffuse radiation predictions because the ground albedo and aerosol-generated 
diffuse irradiance are not included. 

Sunshine hours and cloud-cover data are available easily for many stations: cor¬ 
relations for predicting the insolation with sunshine hours are reviewed. Seven 
models using the Angstrom-Prescott equation to predict the average daily global 
radiation with hours of sunshine have been reviewed. Estimation of monthly aver¬ 
age hourly global radiation was discussed. The 7 models were applied to predict the 
daily global radiation for Hong Kong and comparison with measured data was 
made. Predictions made by Bahel et al. [29], Alnaser [30] and Louche et al. [20] agree 
with the measured data, while the other models provide reasonable predictions. 

Values of global and diffuse radiation for individual hours are essential for many 
research and engineering applications. Twelve decomposition models in the litera¬ 
ture commonly used to predict the diffuse component using the measured hourly 
global data are reviewed. The models were used to predict the hourly diffuse fraction 
for Hong Kong and compared with the model based on Hong Kong measurements 
[22]. Good agreement was found among the model predictions, but lower diffuse 
fractions were predicted by some models [15-17,19,21] for clear-sky conditions. The 
beam-radiation index for Hong Kong was estimated by the Louche et al. [20], Vig¬ 
nola and McDaniels [40] and Maxwell [23] models and compared with the prediction 
of Lam and Li’s model [22]. Again, good agreement was found except that a lower 
beam-radiation index was predicted by Lam and Li’s model [22] under clear-sky 
conditions. 

The clearness index k t and diffuse fraction k ( i for Hong Kong were also deter¬ 
mined with the Iqbal model C [8] for 4 atmospheric conditions (clear, clean, turbid 
and very turbid atmospheres) for 3 cases of ground albedo (corresponding to black 
concrete, uncoloured concrete and white glazed surface). The results are then com¬ 
pared with those predicted by Lam and Li’s model [22]. It was found that Lam and 
Li’s predictions agree with the modeled results for F is about 28 km for k t between 
0.4 and 0.7, and fall for V is between 11 and 28 km for k t between 0.2 and 0.4. The 
parametric model would provide accurate predictions of solar radiation as turbidity 
information is included and good evaluations of thermal environments in buildings 
are therefore possible. If precise atmospheric information, such as optical properties 
of clouds, cloud amount, thickness, position, the number of layers and turbidity, is 
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not available, decomposition models based on measured hourly global radiation 
would be a good choice. 
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